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O ; Abstract 
O 

' Accurate calculation of the gradual insphal motion in an extreme mass-ratio binary system, 

■ in which a compact-object inspirals towards a supermassive black-hole requires calculation of the 

. interaction between the compact-object and the gravitational perturbations that it induces. These 
metric perturbations satisfy linear partial differential equations on a curved background spacetime 

J> ' induced by the supermassive black-hole. At the point particle limit the second-order perturba- 

\o : 

^3 . tions equations have source terms that diverge as r~^, where r is the distance from the particle. 

(N 
O 

\o 
o 



This singular behavior renders the standard retarded solutions of these equations ill-defined. Here 
we resolve this problem and construct well-defined and physically meaningful solutions to these 



I 



O ■ equations. We recently presented an outline of this resolution []|. Here we provide the full de- 

$-H . tails of this analysis. These second-order solutions are important for practical calculations: the 

bX). 

planned gravitational-wave detector LISA requires preparation of waveform templates for the po- 



^ ■ tential gravitational-waves. Construction of templates with desired accuracy for extreme mass- 

ratio binaries requires accurate calculation of the inspiral motion including the interaction with 
the second-order gravitational perturbations. 
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I. INTRODUCTION 



Consider a binary system composed of a small compact-object with mass fi (e.g. a 
neutron star or a stellar-mass black-hole) that inspirals towards a supermassive black-hole 
with a mass M. Such extreme mass-ratio binaries (e.g., M/ fi = 10^) are valuable sources for 
gravitational waves (GW) that could be detected by the planned laser interferometer space 
antenna (LISA) 2]. To detect these binaries and determine their parameters using matched- 
filtering data-analysis techniques one has to prepare gravitational waveform templates for 
the expected GW. An important part in the calculation of the templates is keeping track of 
the GW phase. Successful determination of the binary parameters using matched-filtering 
techniques often requires one to prepare templates with a phase error of less than one-cycle 
over a year of inspiral 0]. Calculating waveform templates to this accuracy is a challenging 
task since a waveform from a year of inspiral may contain 10^ wave-cycles Q]. 

To carry out this calculation one is required to calculate the compact-object inspiral 
trajectory. By virtue of the smallness of the mass ratio fi/M one may use perturbations 
analysis to simplify this calculation. In this analysis the full spacetime metric is represented 
as a sum of a background metric - induced by the supermassive black-hole, and a sequence of 
perturbations - induced by the compact-object. The object's trajectory is also treated with 
perturbations techniques. At leading order of this approximation the object's trajectory 
was found to be a geodesic in the background geometry (see e.g., |5i]). At higher orders, the 
interaction between the object and its own gravitational field gives rise to a gravitational 
self-force that acts on the object. The leading order effect of the self-force originates from 
interaction between the object and its own first-order gravitational perturbations that are 
linear in fi. This first-order self-force (that scales like /i^) induces an acceleration of order 
fi for the object's trajectory. In the case of a vacuum background geometry, formal and 
general expression for this first-order gravitational self-force was derived by Mino, Sasaki, 
and Tanaka [6], and independently by Quinn and Wald [3] using a different method. Later, 
practical methods to calculate this self-force were developed by several authors 
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14, HQ, see also HQ 



;iiiuu. ijaici , 



for a different approach to this problem. The next order 



corrections to the object's trajectory originate from the interaction between the object and 
its own second-order gravitational perturbations (that are quadratic in fi). Higher-order 
corrections will not be considered in this article. 
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We shall now estimate the effect of the gravitational self-force (more accurately the dis- 
sipative part of the gravitational self-force, see below.) on the accumulated phase of the 
emitted GW. This will allow us to determine how many terms should be retained in the 



perturbations expansion (see also [19|, |20|, |21[). For simplicity consider a compact object 



which inspirals between two circular orbits in a strong field region of a Schwarzschild black 
hole. As the object inspirals towards the black hole its orbital frequency slowly changes from 
its value at an initial time. This shift in the orbital frequency is approximated by ujt, where 
uj = and t denotes the elapsed time (from initial time) in Schwarzschild coordinates. Let 
A0 denote the part of the phase shift of the GW (between two fixed times) which is induced 
by the shift in the orbital frequency. Since the GW frequency is proportional to the orbital 
frequency we find that after an inspiral time Aims, the phase shift A0 is approximately 
proportional to Atf^gUJ. Let us find how the quantities in this expression scale with fi. The 
inspiral time Aims scales like AEE~^, here E denotes the particle's energy per unit mass, 
AE is the energy difference between initial and final circular orbits, and E = The 
first-order gravitational self-force produces the leading term in an expansion of E, which we 
denote Ei. Since Ei scales like /i we find that Atins scales like M'^fi^^. Turning now to u, we 
write this quantity as cI; = ^E. At the leading order ^ is independent of yU — it is obtained 
from the equations describing a circular geodesic worldline. At higher orders the conservative 
part of the self-force will produce a correction to this geodesic orbit. Since we focus on the 
contribution to the phase coming from the dissipative part of the self-force (i.e. the part of 
the self-force responsible for a non-zero E.) we ignore the conservative corrections. Denoting 
the leading term in an expansion of u with uji, and recalling that Ei scales like /i we find that 
uji is of 0{fiM~^). Combining the expressions for Atins and ui we find that the first-order 
self- force produces a phase shift A0 of order Atf^^uJi = 0{M/ fi). The second-order self-force 
gives rise to second-order terms in the expansions of E, and uj. These terms are denoted 
here E2 and U2, respectively. Since E2 scales like //^ we find that UJ2 is 0{ij?M~'^). After 
Atins the term uj2 will produce a phase shift of order Atfns^2 = 0[(M//i)°]. Therefore, a 
calculation of A0 to the desired accuracy of order (M/ /i)° (needed for LISA data analysis) 
requires the calculation of the compact-object interaction with its own second-order metric 
perturbations. 

The goal of constructing long waveform templates (e.g. for one year of inspiral) which 
do not deviate by more than one-cycle from the true GW provides practical motivation for 
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the study of the second-order metric perturbations in this article. Moreover, construction of 
second-order metric perturbations allows one to extend the applicability of the perturbation 
analysis to binary systems with smaller M/ /i mass-ratios. This study can also shed light on 
the problem of waveform construction. Suppose that we attempt to construct a waveform 
for an inspiraling compact object (including the correct leading term for A0) by using the 
following procedure. First, we calculate a corrected worldline by including the contributions 
coming from the first-order self-force. Then we substitute this worldline into the expression 
of the source term in the first-order perturbations wave-equation, and finally we construct a 
waveform by solving this equation. Here there is a subtlety, since the first-order gravitational 
self-force is a gauge dependent quantity js^ (e.g. one may set it to zero by an appropriate 
choice of gauge, see below). Therefore, by invoking a first-order gauge transformation we can 
change the path of the corrected worldline. A waveform constructed from this new corrected 
worldline may not encode the correct gauge invariant information. This argument reveals 
that for some gauge choices the above procedure does not provide us with a waveform that 
include the correct leading term for A0. For these gauge choices it is reasonable to expect 
that the correct waveform could be obtained by including the second-order gravitational 
perturbations in the waveform calculation. In this this way all the contributions to the 
waveform which scale like /x^ are being included in the calculation. 

To construct the metric perturbations produced by a compact-object it is useful to con- 
sider the point particle limit - where the dimensions of the compact object approach zero 
(below we give more precise definitions of the compact object and the limiting process that 
we use). In this limit the first-order metric perturbations in the Lorenz gauge satisfy a wave 
equation with a delta function source term (see e.g. j3]). It is well known that (certain 
components of) the retarded solution of this wave equation diverges as as the worldline 
of the object is approached, where r denotes the spatial distance from the object. The con- 
struction of the second-order perturbations is more difficult. In the limit the second-order 
metric perturbation equation away from the compact object take the following schematic 
form 

= v/i(^) v/i(^) & VV/i(^) . (1) 

Here jih^^^ and fi^h^'^^ denote the first-order metric perturbations and the second-order met- 
ric perturbations, respectively. D denotes a certain linear partial differential operator, V 
schematically denotes the covariant derivative with respect to the background metric, and & 
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denotes "and terms of the form...". Since h^^^ diverge as we find that the source term of 
Eq. diverges as r~^. One might naively attempt to construct a standard retarded solu- 
tion to Eq. by imposing Lorenz gauge conditions on h^'^^ and then formally integrating 
the singular source term with the corresponding retarded Green's function. However, the 
resultant integral turns out to be ill-defined, in fact it diverges at every point in spacetime. 
To see this notice that the invariant four- dimensional volume element scales like while 
the source term which is being integrated diverges as r~^. 

In this article we develop a regularization method for the construction of well-defined 
and physically meaningful solutions to Eq. (P). A similar problem in a scalar toy- model 
was recently studied j22]. An outline of the resolution presented in this article was recently 
published yj . Here we provide the complete details of this analysis, including derivations of 

n 

results which were mentioned without derivations in l]. In addition we provide a prescription 



for the construction of Fermi-gauge which was only briefly mentioned in 

For simplicity suppose that the small compact object is a Schwarzschild black-hole. 
We consider the black hole to be "small" compared to the length scales that character- 
ize Riemann curvature tensor of the vacuum background geometry (e.g. a stellar-mass 
Schwarzschild black-hole in a strong field region of a background geometry induced by a 
supermassive Kerr black-hole). Denoting these length scales with {TZi} we express our re- 
striction as yU <^ 7^, where TZ = min{TZi}. The presence of length scales with different orders 
of magnitude allows one to analyze this problem using the method of matched asymptotic 



expansions (see e.g. 



23 



25[). In this method one employs different approximation methods 
to calculate the metric in different overlapping regions of spacetime, where each approxima- 
tion method is adapted to a particular subset of spacetime. Later, one matches the various 
metrics in these overlapping regions, and thereby obtain a complete approximate solution 
to Einstein's field equations. In this article we shall consider the following decomposition 
of spacetime into two overlapping regions. Let r be a meaningful notion of spatial distance, 
we define the internal-zone to lie within a worldtube which surrounds the black-hole and 
extends out to r = rj{TZ) such that rj ^ TZ, and define the external-zone to lie outside 
another worldtube r = r£;(yu), such that fi <^ te- We denote the interior of this inner world- 
tube with S. Since ^ <^ TZ we may choose r^; to be smaller than rj such that there is an 
overlap between the above mentioned regions in r^; < r < r/. We shall refer to this overlap 
region as the buffer-zone (in the buffer- zone r can be of order ^/J7JZ). 
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This article is organized as follows: First in Sec. |n]we discuss the perturbative approxi- 
mation method to Einstein's field equations in the external- zone; in Sec. Ill II we employ this 
approximation and construct the well known first-order metric perturbations; in Sec. llVl we 
discuss the construction of the second-order metric perturbations in the external-zone; in 
Sec. |3we complete the construction of the physical second-order perturbations by matching 
the second-order external-zone solution to a solution in the internal-zone; finally Sec. IVII 
provides conclusions. 



II. APPROXIMATION IN THE EXTERNAL-ZONE 

In the external-zone the spacetime geometry is dominated by the background geometry. 
Therefore, it is convenient to decompose the full spacetime metric (7^"'' into a background 
metric g^^, and perturbations Sg^^i, that are induced by the small black hole, reading 

gif{x) = g^,u{x)+5g^,ix). (2) 

Throughout this paper we use the background metric gf^^ to raise and lower tensor indices 
and to evaluate covariant derivatives. We expand Sg^i, in an asymptotic series, reading 

6g,,{x) = + fi'hf^ix) + 0{fi') . (3) 

Here the perturbations {h^fll} are independent of fi. 

We shall now substitute the asymptotic expansion of (7^"'' into Einstein's field equations 
and obtain linear partial differential equations for the first-order and second-order gravi- 
tational perturbations h^f^} and h^^}, respectively. We assume that full spacetime metric 
satisfies Einstein's field equations in vacuum, reading 

Rif = . (4) 

Here -R^^'' is Ricci tensor of the full spacetime. Substituting decomposition (0) into Ricci 
tensor we obtain the following formal expansion ^ 

Rif = <J + Rii^ M + R^S) M + 0{5g') . (5) 

^ For brevity we shall omit the tensors indices inside the squared brackets of the differential operators. 
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Here the superscripts denote the type of dependence on Sg^y. (0) - no dependence on Sg^j^y , 
(L) - hnear dependence on 6gf^u and (Q) - quadratic dependence on 6g^y. Substituting Eq. 
(0) into Einstein equations ((H) and using expansion (jS)) we obtain 

i?2) = ,x^S ,/ (6) 
Ri'J[h^'^]=0 (7) 

Note that Eq. (0) is an equation for the background metric gf^y. This metric satisfies 
Einstein's field equations in the absence of the small black hole. We can therefore omit the 
restriction x ^ S from Eq. (jH)). 

To further simplify the calculation it is useful to consider the limit /i — of the series 
(j21). Notice that by definition h^^v and h^^u do not depend on /i and therefore the form Eqs. 
(I7I8|) is not affected by this limit. However, the domain of validity of these equations is in 
fact expanded as /i — » 0. At this limit we let r£;(yu) approach zero and Eqs. ()7I8|) become 
valid throughout the entire background spacetime excluding a timelike worldline z{t), where 
r denotes proper time with respect to the background metric. At this limit Eqs. ()7I8|) take 
the form of 

4t)[/.«]=0 ,x^z{t), (9) 

Ri'W^] = -R'^^uW^] ,x^4r)- (10) 

As they stand Eqs. ()9ll0p contain insufficient information about the physical properties 
of the sources that induce the perturbations. To obtain unique solutions to these equations, 
we must provide additional information about these sources. The gravitational perturba- 
tions h^i^u and h^^} are induced by a Schwarzschild black-hole, and therefore their properties 
on the worldline are determined from the physical properties of this source. As we show 
below these properties can be communicated to the external-zone by specifying a set of 
divergent boundary conditions as x ^ -^(t)- Once these divergent boundary conditions are 
specified, a physical solution (defined below) to the perturbation equations (|9ll(jp is uniquely 
determined. In Sec. IVl below we obtain the desired divergent boundary conditions for Eq. 

from the corresponding internal-zone solution. For Eq. (jH)) D'Eath has shown that 
at the limit jj, ^ the retarded perturbations are identical to the retarded first-order 
perturbations that are induced by a unit-mass point particle tracing the same worldhne z{t). 
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Before tackling Eqs. ()9I1U|) we must provide additional information about z{t). Recall 
that as /i — > the worldtube S collapses to the worldline z{t). Roughly speaking one 
may choose S to follow the motion of the black-hole keeping it " centered" at all times with 
respect to 5* in some well defined manner ^. This point of view allows one to identify z{t) as 
a representative "world-line" of the black- hole in the background spacetime. At the leading 
order of approximation the representative worldline of the black-hole was found to be a 
geodesic in the background spacetime (see e.g., 0,0]). Alternatively one may let the black- 
hole drift with respect to the center of S. In this article we find that this alternative point 
of view more suitable for our purposes of studying the second-order perturbations. The 
main reason is that it allows us to choose z{t) to be exactly a geodesic in the background 
spacetime which we denote with zg{t). Setting z{t) = zg{t) guarantees that Eq. Q 
has an exact retarded solution [given by Eq. (fT^ below]. Notice that if we had chosen a 
point of view where z{t) represents the (generically accelerated) motion of the black-hole we 
would have found that Eq. (jH)) does not have any exact solution. This difficulty originates 
form the fact that application of the divergence operator to left hand side of Eq. Q gives 
V^-R^t'' [/i*-^-*] = 0, which restricts the possible sources allowed on the right hand side. As we 
already mentioned, matching with the internal-zone solution implies that one may replace 
the right hand side of Eq. (jUj) with a point particle (delta-function) source at z{t). Here we 
find that the above mentioned restriction on the source implies that z{t) must be a geodesic 
worldline. However, the motion of the black-hole in the background spacetime is generically 
an accelerated motion. This acceleration originates from the gravitational self-force acting 
on the small black-hole. Therefore, had we chosen z{t) to represent the (accelerated) motion 
of the black-hole we would have found that Eq. © does not have any exact solution. Our 
point of view different, since we set z{t) = zg{t) and therefore Eq. © has a well defined 
exact solution. The black-hole acceleration which will gradually shift the black-hole from 
the center of the worldtube S will show up as a term in the boundary conditions for the 
external- zone solution as x ^ zg{t). Notice that since the leading order acceleration scales 
like /i the difference in the boundary conditions induced by this leading-order acceleration 

^ Within the internal zone the fuU spacetime metric is approximated by the metric of a perturbed 
Schwarzschild black hole. Here fixing the center of the black hole amounts to fixing the internal-zone 
dipole perturbations (These dipole perturbations are purely gauge and one can therefore set them to 
zero.) 1^. 
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will affect only n"^v, thus placing the effect of the first-order self-force at the boundary 
conditions for the second-order equation (fTUj). 

With the above choice of worldline the equations for h^j^} and /i^^ now read 

4^,)[/.«] = ,x^za{T), (11) 
i?(t)[/.(^)] = -4'3) [/,(!)] ,x^za{r). (12) 



III. FIRST-ORDER METRIC PERTURBATIONS 

First, we consider the construction of the first-order metric perturbations /t^i\ which 
satisfy Eq. (fTTj). As was previously mentioned, at the limit /i — the perturbations /i^/J are 
identical to the first-order metric perturbations induced by a unit-mass point-particle which 
traces a geodesic zg{t) on the vacuum background metric. We impose the Lorenz-gauge 
conditions on these first-order perturbations, reading 

/i^^)'^^ = , 

where overbar denotes the trace- reversal operator defined by h!]!i} = /i|ti^ — l/2g^i,ha^°' . In 
this gauge the first-order perturbations h^i^} read (see e.g. |25l |) 

/oo 
G^;LpWG{r)K{TV{T)dT. (13) 
-oo 

Here = and (''")] is the gravitational retarded Green's function which is a 

bi-tensor, where the indices refer to z{t), and the indices /x, v refer to x. This Green's 
function satisfies 

UG^\,^,[x\x'] + 2R^%^{x)G^''^,p,[x\x'] = ~A7i-g''^^,{x,x')g%,x')[-g]-^'H\x - x') . 

Here □ = g'^^VpV^ is a differential operator at x, g'^,{x,x') denotes the bi- vector of a 
geodesic parallel transport with respect to the background metric (for the properties of this 
bi-vector see e.g. -^w^ denotes Riemann tensor of the background geometry with 

the sign convention of reference j^], g denotes the determinant of the background metric, 
and 5^(x — x') denotes the four-dimensional (coordinate) Dirac delta- function. Throughout 
this article we use the signature (— , +, +, +) and geometric units units G = c = 1 . 
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IV. SECOND-ORDER METRIC PERTURBATIONS 



We now focus our attention to the construction of the second-order perturbations hliJ 
which satisfy Eq. (jl2|) . Here it will be useful to apply the trace reverse operator to Eq. 
which gives Rl!^^[h^'^^ = —Rl2\h^^^. To simplify the notation we rewrite this equation as 

D^.[h^^^]=S,4h^'^] ,x^zg{t). (14) 

Here we substituted h^^^} = Tij^} — {l/2)g^yhP'''^a into R]!l)[h^'^^] and denoted the resultant 
expression with Z^^^[/i(2)] (notice that RfJ[h^^^] ^ RfJ[h^^^]). Similarfy the source term 
S^v\h'''^^] is defined by S^i,\h^^^] = —R^^\h'^''\ and for abbreviation we shall often simply 
write 5*^1,. The explicit form of these terms is provided in appendix A. For the moment we 
do not impose any second-order gauge conditions and Eq. (fT^ is in a general second-order 
gauge. 

Before constructing a solution to Eq. ()14|) let us first study the singular properties of 
S^v near zg{t). In what follows we shall expand and S^i, in the vicinity of zg{t). 
Throughout this paper such tensor expansions are considered on a family of hypersurfaces 
r = const that are generated by geodesies which are normal to worldline. Each point 
X on such a hypersurface is associated with the same point on the worldline zg{tx)- On 
each of these hypersurfaces the expansions are valid only in a local neighborhood of zg{tx) 
excluding a sphere of arbitrarily small volume which surrounds zg{tx)- Here it should be 
noticed that since the dynamical equations ()11I12|) are not valid on the worldline we do 
not have to keep track of the singularities on the worldline, for example distributions of the 
form 5^{x — zg) which may arise in the tensor expansions below due to application of a 
Laplacian operator on terms of the form 1/r may be completely discarded. Throughout this 
paper (unless we explicitly indicate otherwise) we shall represent the expansions of tensor 
fields using Fermi normal coordinates based on zg{t). These expansions take a particularly 
simple form in these coordinates, and often many of the leading terms are found to be 
identical to the corresponding terms in an expansion over a fiat background spacetime using 
Lorentz coordinates. We shall use the symbol = to denote equality in a particular coordinate 
system. Note that the covariant nature of Eqs. (jllll2p implies that working in a particular 
(background) coordinate system does not reduce the generality of our analysis, since once a 
solution is constructed in one particular coordinate system it can be transformed to any other 
coordinate system by a coordinate transformation. Moreover, our results which provide a 
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prescription for constructing the second-order gravitational perturbations are stated in a 
covariant manner, and therefore can be implemented in any coordinate system. 
Using Eq. (fT^ we expand in the vicinity of zg{t), which gives 

h^;^Xx)=Au,u,r-' + Oir^). (15) 

In the external-zone r denotes the invariant spatial distance along a geodesic connecting 
z(tx) and x, r = \/5^^^bX°^ where are the spatial Fermi coordinates; = (5g is a vector 
field which coincide with four-velocity on the worldline. By substituting Eq. (fT^ into S^i, 
we obtain the following expansion 

S^,{x) = [Au^u, + 7r/^, - Un^n,] r-^ + 0{r-^) . (16) 

Here r]^y denotes Minkowski metric, and we defined Vt"' = x"/r, Vl^ = 0, and Vl^ = g^^Vl^. 
Naively one may try to construct the standard retarded solution to Eq. (fT^ . say by imposing 

(2) 

Lorenz gauge conditions on hiiv and then formally integrating Sij,u with the retarded Green's 
function which gives 

Here there is a problem, examining Eq. reveals that 5*^,^ diverges like in the vicinity 
of the worldline zcij). Recalling that ^y—g{x')d'^x' scales like r^, we find that the integral 
in Eq. ()17|) diverges at every point in spacetime. Furthermore, the next order term in Eq. 
()16|) that diverges like also gives rise to a divergent integral. 

We will now develop a method to obtain well defined solutions for Eq. (fT^ . This method 
is based on consecutive steps, where in each step we reduce the degree of singularity of the 
field equation at hand. Eventually we end up with a field equation for a certain residual 
potential which has a source term which diverges like r~^, this equation has well defined 
retarded solutions. At this point we shall construct a retarded solution to this equation and 
discuss the matching to the internal-zone solution. 

We should mention here another difficulty in calculating the integral in Eq. (|17|). Asymp- 
totically h^^^ has a form of a gravitational wave. Therefore, the leading asymptotic behavior 
of V/i^^) is 0{R~^), where R is the area coordinate - in this paragraph, for simplicity, we 
specialize to a background metric of a Schwarzschild black hole. This implies that for an 
infinitely long world line the source term Sa/3 decays asymptotically as R~^. This term has a 
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static 0{R~^) part which does not vanish after time averaging. An attempt to calculate the 
integral (fT7|) over this static 0{R'^) part produces a divergent integral (even if we resolve 
the difficulty with the singularity near the world line). The regularization of this divergency 
ies outside the scope of this article (Ori has recently suggested a resolution to this problem 



2J|). In what follows we shall assume that such a regularization at infinity has been carried 
out. 



A. r ^ singularity 

First we tackle the strongest singularity in the source term of Eq. (jl4j) . For this 

purpose let us decompose hii,u into two tensor potentials, reading 

/^S = + ' (18) 

where ipfj^u satisfies 

D^4^P] = [Au^u, + I'n^, - Un^n,] r-^ + 0(r-^) . (19) 

Notice that the singular term in Eq. (jl9p is the same as the singular term in the 
expansion of S^y [see Eq. flTH|)]. but the lower order terms of equations (fT^ and (fTUI) are in 
general different. In fact, we do not impose any restrictions on the lower order terms in Eq. 
flTI?|) . Suppose that we construct a solution to Eq. (fT^ . then by subtracting D^^^ip] from 
both sides of Eq. (fT^ we obtain the following equation for 5h\i{, 

D,.,[6h^^^] = S^, - D^,[i,] , X ^ ZG{r) . (20) 

By construction the source term in this equation diverges only like r~^, while the original 
field equation (fT^ has a source term which diverges like r~^. In this sense Eq. (^0]) is 
simpler then Eq. fll4|) . 

We now face the problem of solving Eq. (fT^ . To construct a particular solution we use a 
linear combination of terms which are quadratic in li]^}. Since /i^ti^ diverges like we find 
that by applying the differential operator D^^y to terms which are quadratic in we can 
obtain terms which diverge like r~^. First we construct four independent quadratic tensor 
fields reading 

^% = h^'^%h^S , V^?. = /^^%e^ (21) 

= , 
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These terms can be combined to form a solution to Eq. ()19|) which reads 

= ^ WaV^u + CB<^M-) - ncc^^ + CD<^S.)] . (22) 

Here the constants ca, cb, cc, cd must satisfy 

ca + cb = 1 , cc + cd = 1, (23) 

but are otherwise arbitrary. One may directly substitute Eq. ()22|1 into Eq. and verify 
that ip^i, satisfies this equation. 

Eq. ()22j) can be derived as follows. First notice that the coefficient in front of the 
term in both Eq. (fT^ and Eq. (fT^ does not depend on the curvature of the background 
spacetime. In fact the form of this expression would not change if we would replace the curved 
background spacetime with a flat spacetime. Therefore in deriving -i/^^j^ we may consider a 
simple case of flat background spacetime. In this case, the exact nonlinear solution to 
Einstein's field equations in our problem is simply the Schwarzschild solution. Far from 
the black-hole this Schwarzschild solution may be approximated by an expansion which 
schematically reads 

g':' = V,u + fiH^^J + f^'Hj^J + 0(/i3) . (24) 

_ —————— 

In the appropriate coordinates the second-order term H^J satisfies Eqs. ()14|19|1 in the fiat 
background case. To be consistent with our first-order (Lorenz) gauge conditions we have to 
make sure that the term H^u) satisfies the Lorenz gauge conditions for a fiat spacetime metric. 
As we will immediately show this condition is satisfied if we express the Schwarzschild solu- 
tion in isotropic Cartesian coordinates. The Schwarzschild metric in the isotropic Cartesian 
coordinates takes the form of 




Here = x'^+y^+z^. Expanding this metric in powers of /i/r gives the following expressions 
for the trace-reversed first-order and second-order perturbations 

Hl^J = , ^^i} = - 4^ [2«,S. + rr/,.] . (26) 

Here vP- = (5q is a vector field. Notice that the Lorenz gauge conditions H^^^'^'^^u = are 
satisfied. We may replace h^fj^} with Hju) in the quadratic terms (pT|l and employ Eqs. pHjl 
to express HliJ as a linear combination of these quadratic terms. The coefficients of this 
linear combination are the desired coefficients in Eqs. (j22l23|) . Notice that for the case of a 
fiat background space-time ipfj^^, is an exact solution to Eq. 
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B. r ^ singularity 



We now consider the construction of 6h^^u which satisfies Eq. i^20\i. This equation has a 
source term which diverges hke r~^, and therefore its standard retarded solution diverges. 
Let us examine the terms that give rise to this singularity. It is convenient to express 
the source term of Eq. schematically (and without indices) as 

S - D[ij] = W^^VU^'' & h^^^VVU^'' . (27) 

Notice that here only sum over all the terms diverges as (Since the individual terms 
which di verg e as cancel each other.). Using a decomposition devised by Detweiler and 
Whiting |28[ we decompose h^/^u as follows'^ 

h^=h^' + h^''- (28) 

Here h^^^^ is a certain singular potential which diverges as as r — > 0, and h^f^^^ is a 
certain regular potential which satisfies the following homogeneous wave equation 

+ 2i?V,/.W^ = . (29) 



Decomposition (j28p is particularly useful for expressing the first-order gravitational self 



force, since the general expression of this self-force is completely determined from hi 



see 



28 1 and also Eq. (jH^j) below]. Generically /il/;]^ is a smooth field of 0(r°) on zg{t). 
Expanding hltu and its covariant derivatives in the vicinity of the worldline zg{t) gives (A 



method of constructing these expressions is described in detail in |25l |) 

/i«^ = ^w^M. + 0(ri) (30) 

"^p'^iu' = -^«M«-^P + 0(r°) (31) 

V,V,/iW^ = [3fi,fi, - u,u, - r/J + O(r-i) . (32) 



Notice that the orders r°, r ^ and r ^ are missing from the expansions of hliu , V phiiu and 
Vr,VphitJ , respectively. The absence of these terms can be traced to the vanishing accel- 
eration of zg{t). We now substitute Eq. (j2HI) into Eq. (j2Zj) and examine the various terms 
that give rise to the problematic singularity in the source term S — D[il)]. Expansions 



^ Here we define h^/}}^ and /i^l]"^ to be independent of /z, as opposed to |2 
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imply that the only combination that produces this r ^ singularity is of the form 



To eliminate this problematic singularity we utilize gauge freedom and employ a 
first-order (regular) gauge transformation x"^ —>■ x'^ — fiC,'^, which gives 

^11^ - + ^r^"^'"^ ' + + (33) 

Here does not depend on /i. Notice that we have included the entire gauge transformation 
in the definition of the new regular potential h^f^}^^"'^^\ This is a natural identification since 
the gravitational self-force in the new-gauge is obtained by replacing h^f}}^ with h^f!^}^^"'^^^ in 
the expression of the self-force |29|. We now impose the following gauge conditions 

[hT'^'''%air) ^ [hT + ^z^;^ + ^■^■A.air) = • (34) 
Expanding in the vicinity of the worldline zg{t) gives 

Most beneficially in this new gauge the previously mentioned problematic terms 

source term S — D[ip] diverge only like r ^. This property will 
allow us to construct well defined retarded solution to Eq. (j20|l . Notice that we invoked a 
regular gauge transformation in the sense that it did not change the singular properties of 
hl}u near the worldline, meaning that Eq. is unchanged by this gauge transformation. 
Therefore, the coefficient in front of the term in Eqs. p4ll9|) is not affected by the gauge 
transformation. This implies that even though the numerical values of ip^i, are changed by 
the gauge transformation, the general form of ipfj^t, given by Eqs. ()22l23p is invariant to any 
such regular gauge transformations [e.g., a transformation satisfying Eq. (jBll)]- 



C. Construction of the first-order gauge 

As was previously mentioned the first-order gravitational self-force must be accounted for 
when imposing boundary conditions as x — zg{t) for Eq. (fT^ . To simplify the calculation of 
these boundary conditions we once more use the gauge freedom. The first-order gravitational 
self- force is a gauge dependent quantity js^ , and in fact one can always choose a convenient 
first-order gauge in which the first-order gravitational self-force vanishes (see below). In 
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this gauge the geodesic worldhne zg{t) represents the black-hole's worldline accurately up 
to errors of order /i^. In this case the contributions to the boundary conditions of Eq. 
(|12|) that originate from the black-hole's acceleration due to the first-order self-force simply 
vanish. 

To spell out the desired gauge conditions let us examine the expression for the 0(/i) 
acceleration which is induced by the first-order self-force j^l 

= -fiigf^^ + u^unu''u''i\/,h^^J'' - ^V./i«^) . (35) 

Here all quantities are evaluated on the worldline. Originally this expression was derived for 
h^i^u^ which satisfies the Lorenz gauge conditions. But following the analysis in JsiJ] we find 
that this expression is also valid in any new gauge provided that the gauge transformation 
from Lorenz gauge to this new gauge is sufficiently smooth. To obtain a gauge with a 
vanishing first-order gravitational self-force we should require that a'^ = 0. This requirement 
conforms with many gauge choices. For example it is satisfied if all the first covariant 
derivatives of the regular field in the new gauge vanish. Putting this gauge condition together 
with our previous gauge condition (jH^ yields 

[„a,H,P™.)] ^^^^^ , [„(.)H + j^^^ + j^^J _ _ (38) 
[V^A™"'""""'] ^ [v^Ce'" + {«. + = . (37) 

We shall refer to this new gauge as Fermi gauge. 

To construct Fermi gauge consider contracting Eq. ()37p with u'^. The resultant equation 
states that jg constant along zg{t), which is consistent with Eq. (jHSl)- But more 

importantly it implies that once Eq. (jHUj) is satisfied at an initial point ^^(''0) Eq. (j?fj) will 
guarantee its validity everywhere along zg{t). We now choose an arbitrary gauge vector 
^(o)/i at some initial point zg{to), and construct its first covariant derivatives at this point 
such that Eq. is satisfied. For example we may choose 

To transport and ^fj^-iy along zg{t) we derive transport equations as follows. We treat Eq. 
(|H7|) and the commutation relation 2^^.[^a] = R''f_iua^<i ^ algebraic equations for Ca;/37 

and use the identities of Riemann tensor to obtain the following relation 
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Here all quantities are evaluated on the worldline. One may substitute Eq. (jH^ into Eq. 
(jnZj); and into the above mentioned commutation relation; and thereby verify that these two 
equations are identically satisfied by Eq. (jSHI). We now construct a second-order transport 
equations for .^^(r) by contracting Eq. (j39|) with u°'u'^ which gives 

Solving this equation with the above mentioned initial conditions provides us with the gauge 
vector S,^{t) along the worldline. Similarly we can construct a first-order transport equation 
for ^^.y by contracting Eq. (jH^ with -u". Using the solution ^^(r) of Eq. pUj) together 
with the initial conditions (jnH|) this first-order transport equation can be integrated to give 
iij.;u{,'T)- i^jt-vaij) is then obtained by substituting iptij) into Eq. (jSHl)- Once ^^(t), ^fj,;v{j)i 
and i^-^uai.'^) are obtained, one can use these quantities to construct a local expansion of 
the gauge vector field ^^{x) in a local neighborhood of the worldline. Notice that Eq. (|39p 
implies that Fermi gauge satisfies the Lorenz gauge conditions along the worldline. 

The above construction only provides leading terms in an expansion of in a local 
neighborhood of the worldline. One may continue globally to the entire spacetime. Since 
gauge freedom is associated with non-physical degrees of freedom we are allowed to introduce 
a gauge continuation which depends on arbitrary parameters. Nevertheless, it is sometimes 
helpful to work in a gauge which is manifestly causal, such a gauge may be constructed by 
continuing along future null cones based on zg{t) |31j. In the analysis below we assume 
that such a global gauge continuation has been performed and that the first-order gauge is 
fixed globally. 

D. Particular solution 5h^'^^ 

We shall now construct a particular retarded solution to Eq. (pUj). Here it is useful to 
remove the restriction x ^ zg{t), and continue the source of Eq. pO|l to the world line. 
Clearly a particular solution to Eq. ()2()|1 with the world line included also satisfies the 
original equation (i.e. with the worldline excluded). However, not every continuation of 
the source of Eq. (^01) to the world-line produces an equation which is self-consistent. This 
is easily demonstrated by taking the covariant divergence of both sides of Eq. ()20|) . In 
appendix B we show that V ^ D ^i,\5h'^'^'>] vanishes identically, and furthermore the covariant 
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divergence of the source term of Eq. ()20|) inx ^ -^gI^) vanishes as welL These facts constrain 
the permitted continuations of the source of Eq. ()20|) to the worldhne. If one naively chooses 
a continuation to the worldhne which has a non- vanishing covariant divergence, the resulted 
equation will not be self consistent. Here we choose the simplest possible continuation by 
requiring that no additional singularities are introduced on the worldhne. Meaning that the 
singularities of the continued source term on the worldhne are completely specified by its 
expansion in x ^ zg{t), and no additional singularities (e.g. delta functions) are introduced 
on the worldhne. Eq. (I^Uj) now takes the form 

D,46h^'^]=6S;:,. (41) 

Here SSj^^ = 5^^— Z^^i^f-j/^^], where the superscript F indicates that source terms are evaluated 
in Fermi gauge. Below we show that Eq. (j^T|) is self-consistent by constructing a solution 
to this equation. 

Recall that we have fixed the first-order gauge, but we still have the freedom to invoke a 
purely second-order gauge transformation of the form 

Here the gauge vector ^^-^ is independent of the mass fi. Similar to the first-order case, 
one may choose the gauge vector ^|^2) such that 5h^^u satisfies the Lorenz gauge conditions, 
reading 

5/i(2W;^ = 0. (42) 

Eq. 1)411) now takes the form of 

□<5/^g) + 2RyjU^^ = -2<, . (43) 

We define 5hfj,J to be the retarded solution of Eq. (jlHj) reading 

Sh^i^i^) = ^J G,/^'-*[x|x']5^f,,,(x') x/^^ci^a;' . (44) 

Since the source term of Eq. diverges only like r"^ the integral in Eq. (jl^ has a finite 
contribution originating from the vicinity of zg{t). Notice that even though the expression 
in Eq. ()44j] satisfies Eq. it is not a priori guaranteed that it also satisfies Eq. (jHJ. 

This equation will be satisfied only if the retarded solution (jl^ satisfies the Lorenz gauge 
conditions (021); this is shown in appendix B. 
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E. General second-order solution 



So far we have constructed a particular solution to Eq. (fT^ . reading 

h^;} = + fh^S^ . (45) 

Having found one particular solution does not complete the construction, since we need to 
make sure that the constructed solution satisfies several required physical properties (e.g. 
it has to match the internal-zone solution). To find the desired physical solution we first 
construct the general solution to Eq. (jl4j) . and then impose a set of additional requirements 
on this solution. In this way we obtain a particular solution which is physically meaningful. 

Since Eq. (fT^ is valid for x ^ zg[j\ we find that we can construct a new solution by 
adding to h^l a potential that satisfies a semi-homogeneous equation i.e., a homogeneous 
equation for x ^ zg{j), reading 

I^,.[/^(')''^]=0 ,x^zg[t). (46) 
The general solution to Eq. (fT^ is given by 

- e^"'' + • (47) 

where hii,J is the general solution to Eq. (|4d]). 



F. physical second-order solution 

To find a (particular) physically meaningful solution to Eq. (jl4p we need to impose 
additional requirements on h^^u'^ . These requirements can also be expressed as require- 
ments imposed on the general semi- homogeneous solution h^^u^^ , thus obtaining a particular 
semi-homogeneous solution. For abbreviation we denote this particular semi-homogeneous 
solution with 7^^,. Using Eqs. ()45|47|) the desired physical solution h^^}^ is expressed as 

/^ir = h^;} + 7.. = + Sh^'} + %u . (48) 

We group the additional requirements into four groups (i) gauge conditions (ii) causality 
requirements (iii) global boundary conditions (iv) boundary conditions at the worldline. 
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(i) First we impose gauge conditions. To obtain a simple representation for 7^1, we impose 
the Lorenz gauge conditions on 7^,^. In this gauge Eq. takes the form of 



□7;.. + 2i?y ,7r,p = Q ,x^zg{t). (49) 

Notice that by construction both 5h]iu and 7^jy satisfy the Lorenz gauge conditions. However, 
does not satisfy these conditions, and therefore our particular second-order solution li^ii}^ 
in not in the (second-order) Lorenz gauge. 

(ii) Next we discuss causality. Causality is more easily discussed in terms of an initial 
value formulation of the problem. Therefore, in this paragraph only we consider such an 
initial value formulation. Suppose that we prescribe initial-data for the first-order and 
second-order metric perturbations on some initial spacelike hypersurface So, such that the 
corresponding constraint equations are satisfied on this hypersurface. Here we consider a 
standard extension of the retarded solutions h!]^}, 6h^^u [given by Eq. (fT!?|) . Eq. (jH))] to 
include additional terms which depend on the initial-data. Let x be a point within the 
causal future J~^{Tjq). Then by construction the retarded solution h^^u{x) is unaffected by 
an arbitrary modification of the initial data on outside J^{x) fl Sq. At second-order we 
define 7^1. to be the retarded solution of Eq. (jl^ . Recall that dTijiu is the retarded solution 
of Eq. (|43p . and ip^^, is completely determined from the first-order metric perturbations. Eq. 
(j48j) now implies that the particular solution Tij}}^ is unaffected by an arbitrary modification 
of the initial-data outside J^{x) fl Sq. In this sense the constructed second-order solution 
h^^u^{x) is manifestly causal. 



(iii) In addition we require that the only source for j^iy is the small black-hole. Since 6h 



(2) 

satisfies an inhomogeneous equation ()43|) it may contain waves that are not sourced by the 
worldline. However 7^j, satisfies a semi-homogeneous equation (j49j) . and therefore the waves 
within can originate either from the worldline or from global boundary conditions (e.g. 
prescribed boundary conditions at X^). Since we are interested only in perturbations that 
are induced by the black-hole, we exclude any supplementary perturbations coming from 
these global boundary conditions. 

(iv) We now turn to discuss the boundary conditions as a; — > zg{t). For this purpose 
let us consider once more a small but finite value of fi. In this case Eq. (jf^ is valid only 
for X ^ S. Following D'Eath analysis of first-order metric perturbations 5| we express a 
solution to this equation using a Kirchhoff representation. In this way 7^1, is expressed as a 

20 



certain integral over a surface of a worldtube. Recall that the external-zone lies outside a 
worldtube with radius r£;(/i). Denoting the surface of this worldtube with S^;, we express 
l^iiy as 



Here dTi^/ denotes an outward directed three-surface element on S^;. In the derivation of Eq. 
()50|) we assumed that 7^,^ decays sufficiently fast at spatial infinity. Furthermore, we assumed 
that the retarded Green's function falls sufficiently fast into the past. Consider substituting 
a given expansion of 7^j, (in powers of r) into Eq. (j5(J|) and then taking the limit ^ 0. 
Recall that at this limit te ^ and notice that (iSg scales like r"^. Therefore, only the 
diverging terms (as r ^ 0) in this expansion give rise to a non-vanishing contribution to 
7^,y at the limit /i — 0. We conclude that at the limit, it is sufficient to specify divergent 
boundary conditions to obtain a unique physical solution to Eq. (j49|) . To obtain these 
boundary conditions we examine the divergent behavior of 7^;/(a;) as x ^ zg{t). For this we 
use Eq. (jlH)) together with an analysis of the behavior of -j/^Jj,, and h^^u^ near r = 0. 

First let us consider the divergent behavior of tpj^^. Using Eqs. ()28I30I36|) together with 
Eqs. ()21I22|) we obtain the following expansion for ipj^^^ in Fermi normal coordinates 



Next we examine the behavior of 6hli,J near r = 0. Solving Eq. ()43p iteratively (see Appendix 
C) shows that 6h^^u is bounded as r 0. Finally, we have to determine the divergent 
behavior of in the vicinity of r = 0. This requires an analysis of the internal-zone 
solution, which is discussed in the next section. 

V. APPROXIMATION IN THE INTERNAL-ZONE 

To simplify the discussion we assume that all the length scales characterizing the back- 
ground spacetime {TZi} are of the same order of magnitude TZ. Recall that fi TZ, and 
furthermore in the internal-zone we have r <^ TZ. By virtue of the smallness oi r/TZ and 
fi/TZ we may expand the full spacetime metric in the internal- zone as follows 




p * 



1 

4^2 



(51) 



9tiu = g\ 



Sch 



+ n 



9^J + n-'9i^J + o{n-'). 



(52) 
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Here g^^'^ is the metric of the small Schwarzschild black-hole. Recall that in the buffer-zone 
both expansions Q and are valid, and therefore we can formally expand these two 
equations simultaneously using the fact that in the buffer-zone both yur~^ and rTZ~^ are 
small. Following Thorne and Hartle the various dimensional combinations involved in 
these expressions can be summarized in a table (see Table I). 



TABLE I: Schematic representation of dimensional quantities combinations in expansions of the 
metric in the buffer zone. The top row gives the metric's external-zone expansion and the left 
column gives the metric's internal- zone expansion (|52() . 
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Sch 
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The top row in Table I gives the external-zone expansion, and the left most column in this 
table gives the internal-zone expansion. In the buffer-zone where both expansions are valid 
one can find appropriate coordinates in which these two expansions coincide. Each entry in 
this table schematically represents the combinations of the dimensional quantities (/i,7^, and 
r) obtained from the simultaneous expansions of the top row and left most column. Note 
that this table provides only the powers of the relevant dimensional combinations and does 
not give the exact expressions. 

As before the background metric (7^1/ is described in Fermi normal coordinates based on 
zg{t). The Schwarzschild metric g^^^ is described in the Schwarzschild isotropic coordinates 
(I25|l . rj (schematically) denotes the Minkowski metric, which is the leading order is both the 
expansions of g^^ and g^^^ in the buffer-zone. 

We number the rows from top to bottom starting at the row of g'^^^ (row 0), and number 
the columns from left to right starting from the column of gfj^^, (column 0). Each entry in the 
table is associated with an ordered pair of numbers: (row, column). The divergent behavior 
of /liS'' near the worldline follows from column 2. In this column only terms (0, 2) and (1, 2) 
need to be considered since only these terms have the potential of producing divergent terms 
as r 0. These terms follow from an expansion for g^^^ and TZ'^gjn) that we discuss next. 
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Let us consider first row 0. In the Schwarzschild isotropic coordinates ()25p the entries 
(0, 1) and (0, 2) are nothing but the terms fiHjlJ and jj^H^^y [see Eq.(j2H)], respectively; and 
their exphcit form can be easily obtained from Eq. (|26|). We identify f with r, with this we 
find that the term (0, 1) coincide with the corresponding terms of order r^^ in the expansion 
of . 

We now discuss row 1. Expanding the background metric in the vicinity of the worldline 
gives Qfj^y = rj^i, + 0(r^7?.^^), which implies that the term (1,0) vanishes. This vanishing 
term serves as a boundary condition for the perturbations equations for TZ~^g^^^ as f — oo. 
The term IZ' -i^(^) is obtained by solving the a gravitational vacuum perturbations equation 
in a Schwarzschild background, with these boundary conditions. However, it is well known 
that these 0{TZ^^) perturbations can be eliminated by a gauge choice and mass redefinition. 
Therefore, we may always choose a gauge in which row 1 vanishes identically (see also 
Notice that the vanishing of the (1, 1) term conforms with the fact that the 0{r'^) terms 
are absent from the expansion of in Fermi gauge. Since the (1, 2) term vanishes we 
conclude that the only divergent term in the expansion of h';i' near the worldline, is a term 
which diverges like r~^. The form of this term is provided by H^^u in Eq. ()2fij) . Notice that 
H^^} is identical to the divergent expression of ip^^ given by Eq. (j^H). 



VI. CONCLUSIONS 

— (2')P 

We have found that the divergent terms in the expansion for hiil coincides with the 
divergent terms in the expansion for %l>^^ given by Eq. (j3T| . Rewriting Eq. PHj) as 

and recalling that 5hli,{, is bounded as r 0, we find that at this limit 7^1, is bounded as 
well. Recall that only divergent boundary conditions as r ^ can produce a non-vanishing 
semi-homogeneous retarded solution 7^1,. Since the divergent boundary conditions of Eq. 
vanish, we find that 7^,^ vanishes identically. 
From Eq. ()48j) we finally conclude that the physical second-order gravitational perturba- 
tions in the external-zone are given by 

hf/ = + m ■ (53) 
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Here tpj^^, is given by Eq. ((22)), where the perturbations hj}} are in Fermi gauge, and Sh^u is 
given by Eq. ()44|) . Eq. ()53|) provides a simple covariant prescription for the construction of 
the second-order metric perturbations without any reference to a particular (background) 
coordinate system. 
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APPENDIX A: EXPANSION OF RICCI TENSOR 

The linear and quadratic terms in the expansion of Ricci tensor [see Eq.(jS))] are given by 
(see e.g. ^\) 



-hal3;tj.h°'^.,^ + h°'^ {haf3;fj.u + h^j^-afj — 2h 



We also used the notation S^u[h] = —R^^lh], where on the right hand side h^v is expressed 



(Al) 
(A2) 



with h^y. 



APPENDIX B: Sh^^} SATISFIES LORENZ GAUGE CONDITIONS 



Here we show that the Lorenz gauge conditions are indeed satisfied by the retarded 
solution ()44j) . For this purpose we follow a standard method of deriving differential equations 
for V^6h!]i} . By applying divergence operator to Eq. (j43j) and using a contraction of Bianchi 
identities together with the fact that background geometry is a vacuum spacetime, we obtain 



□(V"5/i2.^) = -2^55"^ 



F 



(Bl) 



We assume that Lorenz gauge conditions (j42j) are satisfied on an initial spacelike hypersurface 
Sr and moreover that [(n°Va)V''(5/i}rJ]s, = where n" is normal to S/. We shall now show 
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that the retarded solution of Eq. ()B1|) vanishes, and therefore satisfies the Lorenz 

gauge conditions as required. 

First, we will show that the source of Eq. ()B1|) vanishes for x ^ zg{t). Consider a metric 
gfj_^ that depends on a small parameter /x, and may be expanded as follows 

= g^.{x) + + Ai'j(x) + 0(1^') . (B2) 

Here gfj.u{x) is not necessarily a solution of Einstein's field equations in vacuum, whereas g^i, 
maintain its definition as a vacuum solution to Einstein's field equations. We now employ 
Bianchi identities reading 

r'^V^G^.^O. (B3) 

Here the contravariant metric satisfies g'^^g^p = 5^, denotes the covariant derivative 
with respect to and G^,^ is Einstein tensor evaluated with this metric. We now employ 
decomposition ()B2|) to formally expand Einstein tensor, and the covariant derivative (for 
rank-2 tensors), giving 



^"/^V„ = V^ + liVt + ^^r^ + o(/i=^) . (B5) 

In these expansions the dependence on /i is only through the explicit powers /i*, and 
F2 denote linear operators (defined on rank-2 tensors), whose explicit form is not required 
here. We now substitute Eqs. ()B4IB5|) into ()B3|) and obtain a perturbative expansion of 
Bianchi identities. The Bianchi identities are valid for any value of /i and therefore the 
individual terms in their expansion in powers of /i vanish identically, yielding the following 
set of identities for an arbitrary tensor fields g^i^} and g'j^} 

V^D,.[^«] = 0, (B6) 
WG^^l + V>i[D^p[g'^%^Q. (B7) 

Here we denoted 

Employing Eqs. (jlllB(ilB7IB8|l we find that the source term of Eq. ()B1|) vanishes for 

X ZG{r). 

25 



Next, we show that the source of Eq. ()B1|) on the worldhne is too weak to produce 
a non- vanishing V^dJi^^l. We shall now estimate the strength of the source Eq. ()B1|) on 
the worldline. Consider a hypersurface of constant time, generated by spacelike geodesies 
which are normal to Zg{t). In this hypersurface we consider a small sphere D{e) of radius 
e, centered at r = 0; and calculate the following three dimensional volume integral over 
^^5S^u inside this sphere, reading 

/ -g';;{zG,x')W'''5S%dV' . (B9) 

JD(e) 

If this integral vanishes then the strength of the source term on the world line is weaker than 
a delta-function source term, and it is too weak to produce a non- vanishing V^Shiiu ■ Recall 
that V^dS^^ vanishes for r 7^ 0. Therefore, we may take the limit e — without changing 
the value of this integral. Explicit expression of V^dS^" reads 

v.<^ = {-gy'^S^s^.'-V^a) + {-ar'^^^^s^'^a) - ^a.^js^'' • (bio) 

We now substitute Eq. (jB10|l into Eq. ()B9|1 . and evaluate this integral using Fermi normal 
coordinates, based on the worldline. Notice that in these coordinates gj^'{zG,x') = 6^' + 
O(r^), dV scales like r^, while the second and third terms in Eq. (jBlOj) scale like and 
r~^, respectively. We therefore find that at the limit e the integral ()B9|) over the 
second and third terms in Eq. ()B10|) vanishes. Substituting the first term in Eq. ()B10|) into 
integral ()B9|1 and using Gauss theorem we find that at for small values of e the integral ()B9|) 
is approximated by 

/ 6S;^^'dEa'. (Bll) 

JdD{€) 

Here dD(e) is the surface of the sphere. Consider an expansion of 55*^"' in powers of r in 
the vicinity of r = 0. Here only terms which scale like have the potential of producing 
a non- vanishing integral at the limit e ^ 0. Using the schematic form (fTTj) one finds that 
only terms of the form h^^'^^VVh^^^^ and Wh^^'^^'Vh^^'^^ produce terms in SSj^""' which scale 
like r~^. We now consider an expansion of h^^^^ and its derivatives, see Eqs. ()30l31l32p 
for the leading terms in these expansions. Examining these equations we see that these 
expansions depends on dimensionless quantities of the form ,VL^ ^rj^i,, higher order terms 
in these expansions also include the Riemann tensor and its derivatives. Dimensional analysis 
implies that the terms in SSj^""' which scale like must be linear in Riemann tensor . The 
integral in Eq. ()B9|) provides us with a vector at r = 0. When integrating over the above 

26 



mentioned local expansions we find that this vector must be composed of Riemann tensor, 
u^, and the background metric. However, in a vacuum background spacetime one can not 
construct form these quantities a non-vanishing vector. Therefore the integral in Eq. ()B9|) 
vanishes. 

In the above calculations we showed that the source of Eq. ()B1|) vanishes for x ^ zg{t), 
and furthermore that a volume integral over this source, which includes the worldline van- 
ishes as well. We therefore find (with the above mentioned initial conditions) that the 
retarded solution to of Eq. ()B1|) vanishes, and therefore 6h^^} satisfies the Lorenz gauge 
conditions. 



APPENDIX C: 6h\:J IS BOUNDED AS r ^ 

We show that the retarded potential 6h]^u given by Eq. (jl^ is bounded as r — ^> 0. Recall 
that the source term in Eq. ()43p diverges like as r — 0. Therefore, the integral in 
Eq. ()44|) converges for x ^ zg{t). In particular 6h^^u is finite on the surface of a worldtube 
which surrounds the worldline at a fixed spatial distance r = r^, where tb << TZ. We now 
consider the solution of Eq. ()43j) within this worldtube. By virtue of the smallness of rTZ~^ 
within this worldtube, Eq. (pSj) can be solved iteratively using the following expansions of 
5h!ii} and g^j.^ 



9,u = rl^.u + 0{n~^). (C2) 

Here again we employ Fermi normal coordinates based on zg{t). Note that the time scale 
in which the source term of Eq. PHj) changes is of 0{TZ) and therefore the leading term 
'^^(o)/ii/ satisfies the following equation 

{5'^'d^d,)5h^^^^^ = -25S^^ . (C3) 

Here x°',x^ denote the spatial Fermi normal coordinates. Eq. ()C3|) is a set of Poisson's 
equations for each tensorial component of ^^(o)^,^- To solve these equations we decompose 
into spherical harmonics centered at r = , reading 

Im 
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Im 

The solution for each spherical harmonics component is given by 

= fj ^,r"p^-{r')dr' + B.T. . (C4) 

Here r> and r< are the larger and smaller terms from the pair {r, r'}, respectively; B.T. 
denotes finite boundary terms coming from the contribution of the surface of the worldtube. 
Expanding in a power series gives 

Eq. ()C4|) implies that in the expansion of only the term fl°°(_2)^~^ gives rise to a 
(logarithmic) divergency in (/)|]^, while all the other terms produce a bounded potential 
at r = 0. Using the schematic expression ()27|1 one finds that the terms in the source of 
Eq. that can possibly contribute to the problematic term fl°°(_2)^~^5 ^ire of the form 

U'^'^^WW^^ and V¥^^'^V¥^^'^ . These terms can be expanded in the vicinity zg{t) using 
expansions ()30I31I32|) evaluated to a higher accuracy (see discussion at the end of Appendix 
B). Dimensional analysis implies that o.°^(_2) has to be proportional to a component of 
Riemann tensor. Since we assumed a vacuum background spacetime the only possible non- 
vanishing candidate for a°°(_2) is R^lau^3u"u^ times a constant. Explicit calculation (using 
MATHEMATICA software) of the constant coefficient of this I = term shows that it 
vanishes, which implies that 6h)^^^^^ is bounded at r = 0. The higher order corrections to 
given by Eq. ()C1|) are smaller than the leading term by at least a factor of r7l~^, 

— C2) 

and are therefore bounded as well. We conclude that 6huu is bounded as r — > 0. 
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of ^'^ on E^- , such that ^'^ decays to zero away from the worldhne. In this way the constructed 
gauge preserve causahty in the following sense: The perturbations on the null-cones for 
T < T~ will remain unchanged if one modifies the worldline for r > r~. Such a modification 
of the worldline is possible by introducing additional GW that interact with the worldline for 
r > T~. 
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